On the Anomalous Scaling Exponents in Nonlinear Models of Turbulence 
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We propose a new approach to the old-standing problem of the anomaly of the scaling exponents 
of nonlinear models of turbulence. We achieve this by constructing, for any given nonlinear model, 
a linear model of passive advection of an auxiliary field whose anomalous scaling exponents are the 
same as the scaling exponents of the nonlinear problem. The statistics of the auxiliary linear model 
are dominated by 'Statistically Preserved Structures' which are associated with exact conservation 
laws. The latter can be used for example to determine the value of the anomalous scaling exponent 
of the second order structure function. The approach is equally applicable to shell models and to 
the Navier-Stokes equations. 



The calculation of the scaling exponents of structure 
functions of nonlinear turbulent velocity fields remains 
one of the major open problems of statistical physics 
yj. Dimensional considerations appear to fail to pro- 
vide the measured exponents, and present theory cannot 
even specify the mechanism for the so called "anomaly" , 
i.e. the deviation of the scaling exponents from their 
dimensional estimates. Theoretical attempts to calcu- 
late the exponents were mainly based on perturbative 
expansions or on closures of the infinite correlation 
function hierarchy Q . The aim of this Letter is to pro- 
pose a new idea to ascertain the anomaly of the scaling 
exponents in turbulence. In addition, we exhibit an alter- 
native way to determine the anomalous scaling exponent 
of the second order structure function. The proposed 
approach is equally applicable to Navier-Stokes turbu- 
lence and to simplified models of turbulence, like non- 
linear shell models. The only distinction is in the ease 
of numerical demonstration; for shell models we present 
adequate numerical confirmation of the proposed theory. 
For Navier-Stokes turbulence we present calculations at 
a resolution of 128 3 . 

The central idea is to construct a linear model whose 
scaling exponents are the same as those of the nonlinear 
problem. In this linear problem the exponents are uni- 
versal to the forcing, and we understand the mechanism 
for the anomaly of the scaling exponents; we use this to 
show that also the nonlinear problem must have anoma- 
lous exponents. We exemplify the idea first in the context 
of the Navier-Stokes equations. Consider a model for two 
coupled vector fields u and w, 



du 

~dt 
dw 

~dt ' 



u ■ Vw + Xw ■ Vit = — Vp 



u ■ Vio + Xw ■ 'Vw = —'Vp- 



vV 2 u + f, 
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Here v is the kinematic viscosity, p and p are pressure 
fields imposing V • u = V • w = and / and / are two 
uncorrelated Gaussian random forcing. Finally A is a real 
number. Here and below we assume that the scaling ex- 



ponents are universal to the forcing. We want to demon- 
strate their anomaly and to find their numerical values. 
For A = Eq. Q reduces to the Navier-Stokes equations 
for it, whereas Eq. (J2Jl becomes a linear equation for w, 
passively advectcd by u. This linear problem was referred 
to before as a "passive vector with pressure" It is 

known to exhibit anomalous scaling with exponents that 
are universal to the forcing. In addition , on e understand 
the mechanism for the anomaly 0, 0, fToL ITlj : the lin- 
ear model possesses "Statistically Preserved Structures" 
(SPS) which are evident in the decaying problem. These 
are left eigenfunctions of eigenvalue 1 of the linear propa- 
gator for each order (decaying) correlation function, and 
see below for more detail. 

Evidently, for any finite value of < A < oo the scaling 
exponents of the two fields u and w must be the same, 
due to the symmetry Xw u and the assumed univer- 
sality to the forcing. Consider the two composite fields 
u+ = u + Xw and «_ = u — Xw. Choose the forcing 
terms in such that ((/ + Xf )(f - A/)) = 0. Then 
u + satisfies precisely the Navier Stokes equations (with 
forcing / + A/) while u_ is a passive vector advected 
by u+ and forced by / — A/. Our main proposition 
is that the scaling exponents of the Navier Stokes field 
u + and of the passive vector it_ are the same. If true, 
we can study the anomalous scaling of the Navier-Stokes 
problem by using the successful tools and the concepts 
employed to understand the anomalous scaling for pas- 
sive fields (scalar or vector) . Note that the identity of the 
scaling exponents of u + and u_ is equivalent to saying 
that in Eqs. Q and © the scaling exponents of u and 
w are the same for all A including A = 0. 

To verify our proposition we present in Fig. Q the 
scaling properties of u + and m_ for A = 1 obtained by 
direct numerical simulations of Eqs. and J2J- The 
pseudospcctral code used for the simulations has a res- 
olution of 128 3 . In Fig. (JTJ we show the ESS plot of 
the sixth order structure functions for both w_ and u + . 
One observes convincing scaling behavior with the same 
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FIG. 1: Log- log plot of the sixth order structure functions 
of the fields u+ and U- (circles and squares respectively), for 
A = 1, as a function of the third order structure functions. 
The dashed line corresponds to the best fit in the scaling 
region with slopes 1.77. Lower insert: Sa(r)/ Ss(r) 1,77 . Upper 
insert: z p = ( p /(s— p/3 computed for the structures functions 
of w+ (line) and it_ (circles). 



exponent for both field. In the upper insert we show the 
anomalous exponents z p = Cp/C3 — P/3 for the field u + 
(line) and it_ (circles) computed up to order 8: the agree- 
ment is excellent. Thus, we can indeed propose that the 
Navier-Stokes field u + has the same scaling exponents 
as the passive vector field u_. For additional strong ev- 
idence we turn to shell models 0, 0, Q . 

To reach a deeper understanding of the relation be- 
tween the nonlinear and the linear models, and to clearly 
present the role of the statistically Preserved Structures, 
we consider the Sabra shell model which, like other shell 
models of turbulence, is a truncated description of the dy- 
namics of Fourier modes, preserving some of the structure 
and conservation laws of the Navier-Stokes equations: 

(— + vk^)u n = i(k n+ iu* n+1 u n+2 - <Sfc„U*-l w n+l 
at 

+ (1 - 5)fc„_lli n _lU„_ 2 ) + fn ■ (2) 

Here u n are the velocity modes restricted to 'wavevec- 
tors' k n = fco^t™ with ko determined by the inverse outer 
scale of turbulence. The model contains one additional 
parameter, 5, and it conserves two quadratic invariants 
(when the force and the dissipation term are absent) for 
all values of S. The first is the total energy J2 n \ u n\ 2 and 
the second is J2 n (~ l)™^n \ u n\ 2 , where a = log^l — 5). 
In this Letter we consider values of the parameters such 
that < S < 1; in this region of parameters the second 
invariant contributes only with sub-leading exponents to 
the structure functions [2 113 ■ The scaling exponents 
characterize the structure functions: 

S 2 {k n ) = (u n u n ) ~ k~ G , (3) 
S 3 {k n ) = 3(u„_iu„u* +1 ) ~ fc~ C3 , (4) 

etc. for higher order S p (k n ) ~ fc„ p • 

The values of the scaling exponents were determined 
accurately by direct numerical simulations. Besides £3 



which is exactly unity all the other exponents £ p 
are anomalous, differing from p/3. It was established 
numerically that the scaling exponents are universal, 
i.e. they are independent of the forcing /„ as long 
as the latter is restricted to small n Despite of 

the much simpler structure of Shell Models in compar- 
ison with Navier-Stokes equations no systematic break- 
throughs on the analytical calculation of scaling proper- 
ties has been achieved. Previous attempts being manly 
based on stochastic closures 0, . 

Consider then a passive advected field which in the dis- 
crete shell space has the complex amplitudes w n . The dy- 
namical equations for this field are linear and constructed 
under the following requirements: (i) the structure of 
the equations is obtained by linearizing the nonlinear 
problem and retaining only such terms that conserve the 
energy; (ii) the resulting equation is identical with the 
sabra model when w n = m„; (iii) the energy is the only 
quadratic invariant for the passive field in the absence of 
forcing and dissipation. These requirements lead to the 
following linear model: 



= 3 $4", w) - vk n w n + f n 
where the advection term is defined as 



(5) 



$ n {u,w) = k n+ i[(l + 6)u n+2 w* +1 + (2 - 6)u* l+1 w n+2 ] 
+k n [(l - 25)u* n _ 1 w n+ i - (1 + <5)m„ + iw*_ 1 ] 
+fc n _i[(2 - 5)u„-iw n -2 + (1 - 2S)u n -2W n -i] (6) 

Observe that when w n = u n this model reproduces the 
Sabra model, and also that the total energy is conserved 
because Y^ n 3[$n(", w ) w n\ = 0. The second quadratic 
invariant is not conserved by the linear model. Finally, 
both models have the same 'phase symmetry' in the 
sense that the phase transformations u n —> u n exp (i(f> n ) 
and w n — y w n exp (i6 n ) leave the equations invariant iff 
4> n -i + 4> n = B n -\ + n = 8 n +i- This identical 

phase relationship guarantees that the non- vanishing cor- 
relation functions of both models have precisely the same 
forms. Thus for example the only second and third cor- 
relation functions in both models are those written ex- 
plicitly in Eqs. © and (gj. 

As already remarked, the anomalous scaling of w n can 
be investigated in terms of the SPS F° r exam- 

ple for the second order correlation function denote the 
propagator Pj^ni(t\to)] this operator propagates any ini- 
tial condition {w n w^)(to) (with average over initial con- 
ditions, independent of the realizations of the advecting 
field u n ) to the decaying correlation function (with aver- 
age over realizations of the advecting field u n ) 



K<)(t) - P™,(t\to){w n ,w*,)(t ) 



(7) 



(2) 

The second order SPS, Z„ , is the left cigenfunction with 
eigenvalue 1, 



,(2) 



(8) 



3 




iog(tg 

FIG. 2: The sixth order structure function of the field w n 
in Eqs. JTOJ for A = 1CT 1 , 1(T 3 and 1CT 5 , together with the 
sixth order structure function for the Sabra model and 
for the linear model ©, respectively. The structure function 
of the field u n for A > are not shown since they are indis- 
tinguishable from those of the w n . Inset: log- log plot of the 
fourth-order correlation function ^2,2(^,^7) vs. k n calcu- 
lated for the linear field (+) and for the nonlinear field (solid 
line) at A = 0. 
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FIG. 3: With the symbols (+) the constants 7^ 2 ' (bottom) and 
1^ (top) constructed by projecting the decaying structure 
function of the linear model on the forced structure function 
of the nonlinear model. To emphasize the importance of 
using the correct SPS, we also show the result for i' 4 ' using 
the dimensional Kolmogorov prediction for Z 4 (small dots) 
and Z 4 — 1 (solid line) 



Note that Zn is time independent even though the oper- 
ator ' , (i|to) is time dependent. Each order correlation 

function is associated with another propagator (t\to) 
and each of those has an SPS, i.e. a left eigenfunction 
of eigenvalue 1. These non-decaying eigenfunctions 
scale with k n , ~ k n ^ p , and the values of the ex- 

ponents £ p are anomalous. Finally, one can show that 
these SPS are also the leading scaling contributions to 
the structure functions of the forced problem (JjJJ 0, ^| . 
Thus the scaling exponents of the linear problem 
are independent of the forcing /„, since they are 
determined by the SPS of the decaying problem. 

Using equations JSJ, we can now write the (Sabra) shell 
model version of Eqs. 0,10): 

~dT = + y ^n{w,u) - vk 2 n u n + f n , (9) 

~dY ^ ^® n ( u,w } + y*n( tu ' u; ) ~ v k 2 n w n + /„(10) 

For A = we recover the equations for the nonlinear 
and a linear models, Eqs. J2Jl and (J3J. At this point 
we present strong evidence that the scaling exponents of 
either field exhibits no jump in the limit A — > 0. Accord- 
ingly, the scaling exponents of either field can be obtained 
from the SPS of the linear problem. 

Eqs. 10 and (|10f> were solved numerically, choos- 
ing /„ a constant complex number limited to n = 0,1, 
and /„ a random force with zero mean, operating on 
the same shells. We chose v = 10~ 8 , S = 0.6 and 
A = HT 1 , 10~ 3 , 10~ 4 , 10~ 5 , 0. In Fig. Hwc show, for ex- 
ample, results for the sixth order objects (\u n _iu n u* l+1 \ 2 ) 
and {\w n -iw n w* +1 \ 2 ) . Plotted are double-logarithmic 
plots of these object as a function of k n . We see that the 
exponents of the linear and nonlinear model at A = arc 



the same and they do coincide with the exponents of the 
two coupled models @, (fTUfl for A > 0. The same results 
have been checked numerically up to exponents of order 
10. Hence, the limit A — > is regular. 

We stress at this point that the two problems do not 
share exactly the same statistics; the linear problem, be- 
ing symmetric in w n — > — w n has an even probability 
distribution function (pdf) and thus zero prefactors for 
all the odd structure functions. The statement is only 
about the identity of the scaling exponents, neither the 
trajectory in phase space nor the the pdf. In the in- 
set of Fig. |21 we also demonstrate that the linear and 
the nonlinear problems share the same scaling proper- 
ties for correlations that depend on more than one shell. 
The data pertain to F Ptq (k n ,k m ) = (|it n | p |it m |' 3 ), with 
p = 2, q = 2 for both models. Finally, we comment that 
the limit A — > can be considered mathematically for 
the shell model equations and (|10|) , to prove that it is 
not singular. Such a proof is however beyond the scope 
of this Letter, and will be presented elsewhere. 

The greatest asset of the present approach is that we 
can now forge a connection between the SPS of the linear 
model and the forced correlation function of the nonlin- 
ear problem. This underlines the anomaly of the scaling 
properties of the latter model, and allows us to deter- 
mine £2- We start with the second order quantities. We 
can project a generic second order decaying correlation 
function of the linear model onto the second order SPS, 
thus creating a statistically conserved quantity: 

1^ ^ ^z( 2 )KO(i) = ^z( 2 )p^(t|f )K,<,)(t ) 

= XX'Vn'<'}(*0) • (11) 
n' 
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Where the average is over different initial conditions for 
the linear fields and different realization of the advect- 
ing velocity field. To show that the forced second order 
correlation of the nonlinear field is dominated by Z^ 2 \ 
we use this forced correlation function instead of in 
Eq. IHJ. The test is whether remains constant on 
a time window which increases with Reynolds. This is 
shown in Fig. [31 The success of this test demonstrates 
that (i) there exists a SPS for the linear problem; (ii) the 
SPS is well represented by the forced nonlinear second 
order correlation functions. This is a direct demonstra- 
tion that the correlation function of the nonlinear model 
scales with the same anomalous exponent as Z^ 2 \ An 
even more stringent test can be made using SPS of or- 
ders large than 2, where also correlations between differ- 
ent shells are relevant for the decaying properties [IcHlH . 
For example 1^ is given by the weighted sum of three 
contributions: 

I {A) = En,™ Z^{\ Wn \ 2 \w m \ 2 )(t)+ (12) 

£„ [Z^{w n w 2 n+1 w* n+3 )(t) + c.c.] + 

E„ i Z n A) (w n W n+1 W n+3 W* n+4 )(t) + CO.] , 

where all the terms allowed by the phase symmetry were 
employed. In Fig. we show results for 1^ where 
again we swapped the SPS of the linear problem for the 
measured forced correlations of the nonlinear problem: 

Zn ra — > (|M n | 2 |w m | 2 ) and the corresponding expressions 
for Zn' 4 ^ and Zn ■ We thus conclude that the scaling 
exponents of a given nonlinear shell model can be un- 
derstood from the SPS of an appropriately constructed 
linear problem. To make this point crystal clear, we have 
used in fact the forced structure functions of the nonlin- 
ear model as approximants for Z^ 2 \Z^ in the calcula- 
tion of /( 2 ) and l' 4 ' shown in Fig. [3J The constancy 



of both demonstrates that the forced correlation func- 
tion of the nonlinear model arc very well approximated 
by the SPS of the linear model. This demonstration can 
be repeated with higher order correlation functions with 
the same (or better) degree of success. Finally, the exis- 
tence of a conserved quantity I*- 2 -* can be used to calculate 
£2 = C2 ■ Starting from a given arbitrary initial condition 
(say a 5-function on one shell) and computing Eq. i|ll|> 
with many realizations of the advecting velocity field, one 
finds that there exists a sharply defined £2, ~ k~&, 
for which 1^ is indeed constant. The same approach 
can be used to determine £3 but we know that C3 = 1- 
Unfortunately, this simple approach cannot be used for 
higher order exponents, because the corresponding SPS 
depend on more than one k n , and cannot be represented 
as a simple power law. 

In conclusion the anomalous scaling of nonlinear modes 
of turbulence, be then the Navier-Stokes equations or 
shell models, are fixed by the eigen functions of the in- 
ertial operator, which are precisely the SPS of the linear 
problem. Thus, although the concept of eigenfunctions 
cannot be applied directly in nonlinear problems, we are 
able to argue that the mechanism leading to anomalous 
scaling in Navier-Stokes equations and other nonlinear 
models is identical to the one recently discovered for the 
case of passively advected fields. This conclusion may 
open the way to a deeper understand of intermittency 
in turbulent flows and to a direct computation of the 
anomalous exponents. 
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